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Summary 


This  report  reviews  in  detail  our  progress  on  understanding  the 
short-term  statistics  of  microthermal  turbulence  fluctuations  and  optical/ 
infrared  scintillations,  and  their  interrelationship.  A complete  analytical 
understanding  of  microthermal  and  scintillation  data-spread  and  confidence 
intervals  has  been  achieved,  including  the  effects  of  microthermal-envelope 
probability  distributes,  integral  scales,  varying  averaging  times,  and 
turbulence  intermittency . A practical  definition  and  model  of  intermittent 
turbulence  is  given,  and  it  is  shown  that  the  usual  mathematical  techniques 
may  be  applied.  The  fundamental  determining  factor  in  data  spread  is 
shown  to  be  the  two-point  correlation  of  the  microthermal  envelope  fluctua- 
tions along  the  propagation  path,  with  applications  to  the  short-term 
performance  of  imaging  as  well  as  illumination  systems.  Supporting  data 

are  given. 

As  an  addendum  to  the  preceding  report  on  this  program,  which  deals 
with  turbulence  effects  on  finite  aperture  laser  illumination  systems, 
the  effective  aperture  size  of  an  arbitrarily  truncated  gaussian  beam  is 
derived.  In  addition,  it  is  shown  from  the  Huygens-Fresnel  formulation 
that  the  normalized  variance  of  target  irradiance  for  a large  transmitter 
or  small  coherence  scale  (pq)  is  unity,  in  agreement  with  recent  data  and 
physical  reasonings. 
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I.  INTRODUCTION 


This  report  describes  progress  on  two  projects,  one  in  detail.  In 
the  next  section  (II),  we  describe  at  length  the  theoretical  and  prelimi- 
nary experimental  results  of  our  efforts  on  the  short-term  statistics  of 
turbulence  and  optical/infrared  propagation.  In  the  following  section  (III), 
some  recent  developments  in  the  work  on  finite,  wander-tracking  laser  illumi- 
nation systems  are  briefly  discussed. 


.SHORT-TERM  TURBULENCE  AND  PROPAGATION  STATISTICS 

Significant  conceptual,  analytical,  and  experimental  progress  has 

been  made  on  the  description  of  the  short-term  statistics  of  turbulence 

and  scintillation  phenomena.  In  the  present  section,  this  progress  will 

be  reviewed  in  a largely  self-contained  manner,  such  that  little  or  no 

familiarity  with  earlier  reports1  will  be  assumed  on  the  part  of  the 
reader. 


The  basic  topics  of  interest  here  are  the  description  of  basic  micro- 
thermal  fluctuations,  the  short-term  statistics  of  appropriate  averaged 
quantities,  short-term  measures  of  propagation  effects  and  their  relation- 
ship to  the  micro thermal  quantities,  the  influence  of  the  choice  of  aver- 
aging times,  intermittency  of  turbulence,  and  related  items.  These  aspects 
will  be  discussed  in  the  following  subsections. 

A'  j-tatlst:ics  of  the  Basic  Microthermal  Fluctuations 
The  basic  microthermal  fluctuations  associated  with  the  turbulence 
effects  of  interest  are  known  to  be  characterized  by  probability  distri- 
butions which  differ  markedly  from  normality.  In  particular,  if  the 
temperature  difference  fluctuations  between  two  probes  (AT  ) are  ob- 
served, where  the  probe  separation  is  within  the  inertial  subrange,  the 
distribution  is  symmetric  with  a ‘flatness  factor  which  may  be  as  high 
as  40,  compared  to  a value  of  3 for  the  normal  distribution.  A typical 


to^he^uare  of^he  ^ °f  the  f°Urth  central 
fMr-H  T ? f h SeCOnd  central  moment.  The  skewness  is  the  squared 
third  central  moment  divided  by  the  cubed  second  central  moment.  q 


1.  J.  R.  Kerr,  "Propagation  of  Multiwavelength  Laser  Radiation  Through 
Atmospheric  Turbulence",  RADC-TR-73-322 , August,  1973.  (AD769  792) 
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record  of  such  fluctuations  is  shown  in  Fig.  1,  and  the  probability 

distribution  in  Fig.  2,  exhibiting  a flatness  factor  in  this  case  of  fi/ 

Physically-,  this  exaggerated  tail  on  the  probability  distribution  relates 

to  occasional  fluctuations  of  extreme  magnitude  ("small-scale  intermittency") , 

2 

and  is  fundamental  to  turbulence  phenomena. 

Although  the  above  distribution  is  basically  amorphous,  we  are  more 
interested  in  an  appropriate  measure  of  the  strength  of  turbulence,  and 
therefore  in  the  statistics  of  an  appropriate  average  or  smoothed  quantity. 
This  quantity  is  the  (short-term)  average  of  the  square  of  the  temperature- 
difference  fluctuations,  which  we  obtain  here  through  digital  low-pass 
filtering  techniques. 

3 

Theoretical  treatments  have  predicted  that  non-negative,  differential 
(or  finite-difference)  turbulence  variables  will  exhibit  log  normal  behavior 
wnen  averaged  over  a spatial  or  temporal  scale  corresponding  to  the  inertial 
subrange.  Probability  distributions  of  the  (log)  squared  temperature  dif- 
ference fluctuations  are  shown  in  Fig.  3 for  several  smoothing  or  averaging 

2 

scales,  and  log  normality  is  verified  except  at  small  values  of  (AT^)  • 

This  departure  at  small  values  is  not  due  to  instrumental  noise,  and  has 
been  observed  in  nearly  all  experiments  on  turbulent  temperature  and 
velocity  fluctuations.  The  apparent  effect  of  the  inner  scale  of  turbulence 
is  also  seen  in  Fig.  3. 

As  the  averaging  time  or  scale  size  is  increased,  one  expects,  based 
on  the  central  limit  theorem,  to  observe  the  approach  of  the  distribution 

( 4 

to  normality.  This  approach  is  conveniently  displayed  on  a "beta  diagram", 
which  is  a plot  of  the  flatness  factor  (3~)  vs.  the  skewness  (3.).  The  beta 

2 ^ i 

diagram  for  (AT12)  over  a range  of  averaging  times  is  shown  in  Fig.  4.  Also 
shown  is  the  point  (32  = 3,  3^  = 0)  corresponding  to  normality,  and  the  line 
representing  the  family  of  log  normal  distributions.  The  displacements  of  the 
experimental  points  from  the  log  normal  line,  especially  at  small  averaging 
times,  is  due  in  part  to  the  inability  to  measure  the  higher  moments  of  a 


2.  A.M.Yaglom,  Soviet  Physics-Doklady  LI,  20-29,  July  1966. 

3-  A.S.Gurvich  and  A.M.Yaglom,  Phys.  of  Fields  Sjppl. , Boundary  Layers  and 
Turbulence,  10,  Part  II,  S59-S65  (1967). 

4.  G.J.  Hahn  and  S.S. Shapiro,  StatisticaJ  Models  in  Engineering,  Wiley  & 
Sons,  197  (1967). 
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large-flatness-factor  process  with  suitable  accuracy.^ 

In  order  to  relate  these  short-term  measures  of  mean-square  temperature 
fluctuations  to  the  "strength  of  turbulence"  as  utilized  in  propagation 
theory,  a minimum  averaging  time  is  in  fact  required.  Roughly  speaking, 
this  averaging  time  (spatial  scale)  must  correspond  to  the  lowest  frequencies 
(largest  scales)  in  the  turbulence  fluctuations  which  are  of  interest  in 
the  propagation  phenomena.  For  the  case  of  scintillations  (but  not  neces- 
sarily for  phase  effects),  ve  may  designate  this  critical  scale  as  the 
outer  scale  of  the  inertial  subrange  (Lq).  If  we  average  over  a smaller 
scale,  we  do  not  obtain  meaningful  measures;  if  we  average  over  much  larger 
scales,  we  lose  details  of  interest  in  the  short-term  problem.  The  mean- 
square  temperature-difference  fluctuations,  for  averaging  times  (t)  >_  (L  /v) , 

where  v is  the  wind  speed,  will  be  referred  to  in  this  report  as  the 
2 2 

"short-term  C " or  C . It  is  apparent  from  Fig.  4 that  the  statistics 
n nT  2 

of  a series  of  measurements  of  with  a given  t will  be  log-normal-like, 
and  will  approach  normality  only  If  i is  sufficiently  long,  which  is  not 
the  interesting  case. 

We  note  as  a side  point  that  the  integral  of  any  stationary,  ergodic 
process  is  known  to  be  asymptotically  gaussian.  This  is  true  regardless 
of  the  probability  distribution  of  the  process,  and  may  be  taken  as  a 
generalization  of  the  Central  Limit  Theorem.5  6 

Further  statistical  properties  of  the  mean-square  fluctuations  will 
be  discussed  in  later  sections. 

B.  Integral  Scales  for  Averaged  Quantities 

As  is  apparent  from  the  preceding  discussion,  we  will  be  generally 
concerned  with  random  processes  which  are  defined  as  short-term  averaged 
or  smoothed  versions  of  more  fundamental  random  processes.  An  important 
parameter  for  any  random  process  is  its  "integral  scale",  which  is  a 
meajure  of  the  decorrelation  time  for  the  process.  In  this  section,  we 
will  show  that  the  integral  scale  for  a smoothed  process,  where  the 
averaging  time  is  greater  than  the  characteristic  time  scales  of  the  original 
process,  approaches  (1/2)  the  averaging  time  itself.  That  this  is  so  can 


5.  H.  Tennekes,  j.  C.  Wyngaard,  J.  Fluid  Mech.  55,  93-103  (1972). 

6.  J.  L.  Lumley,  "Stochastic  Tools  in  Turbulence",  Acad.  Pr.,  1970. 
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be  seen  intuitively  by  noting  that  smoothing  corresponds  to  low-pass- 
filtering, and  if  the  filter  cut-off  frequency  is  below  the  characteristic 
frequencies  of  the  original  process,  then  the  lowest  significant  frequency 
(or  reciprocal  integral  scale)  in  the  smoothed  process  will  be  the  filter 
cut-off  frequency  (Fig.  5). 

To  show  this  formally,  we  generalize  the  development  of  Ref.  7.  We 
define  p(t')  as  the  normalized  autocorrelation  or  autocovarlnnce  function 
of  the  basic  (stationary,  ergodic)  random  process  f(t),  and  write  the 
integral  scale  of  f as 

00 

lc  = [ p(t')dt'  c 


We  then  define  the  smoothed  random  process  gT(t)  as 


f (t+t' ' )dt*  1 


where  the  averaging  time  is  t.  For  simplicity,  we  assume  that  f,g  have 

2 

zero  mean,  and  defining  a as  the  variance  of  f,  we  write  the  auto- 
covariance of  g^  as 

1 / fxl 2 rT/2  \ 

<gT  (t+f)  gT  (t)>  --j/j  f(t+t”+t')dt”  f f(t+t'")dt"'\ 

T '-t/2  -t/2  ' 

2 t/2  t/2 

f f dt,,dt,,,p(t,+t,,-t,,,) 

T -t/2  -t/2 

2 t/2  T/2-t 

= £2  f dt  f d t ' ' p (t  '+t 1 ' ) 

T -t/2  -T/2-t 


We  then  integrate  by  parts: 
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(gT(t+t’)  gT(t)) 


p(t,,+t,)dt" 


t/2 


-t/2 


5^  Vt/2  J p (t ' '+t 1 )dt  ’ ' + t/2  J p(t,,+t,)dt' 


-r 


J (t,,-T/2)p(t’-t,,)dt’ 
o 


J (T/2-t,,)p(t,+t”)dt” 


^2  j T^2  / p(t’-t)dt  + t/2  / p(t'+t)dt 


o 

T 


+ J (T/2-t)p(t’-t)dt  + y (T/2-t)p(f+t)dt 


(T-t)[p(t’-t)  + P(t'+t)]dt 


T 


/ 


(l-t/T)  [ P(t'+t)  + P(t’-t)] 


dt 
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As  c.’  special  case  we  note  that 
2 T 

^gT2}  “ ~j  (l-t/i)p(t)dt 

o 

Also,  we  integrate  the  autocovariance  of  to  obtain: 

°°  2 T 
J <gT(t+f)gT(t)>dt’  - ^_If  J (l-t/x)dt 

o O 


= o2I, 


The  integral  scale  of  g^  is  thus  given  by  (5)  and  (4)  as 


o2I, 


<0 


r£  T 


2 J (l-t/x)p (t)dt 
o 


Finally,  if  If  <<  T,  then  the  integral  is  equal  to  If,  and 


lim  I = t/2 
8t 


(4) 


(5) 


(6) 


(7) 


which  is  the  desired  result. 

Typical  temperature-difference  fluctuations  with  various  smoothing 
times  (t)  are  shown  in  Fig.  6,  and  the  autocorrelation  functions  in  Fig.  7. 
In  Fig.  8,  the  computer  integral  scales  are  shown  as  a function  of  t,  and 
it  is  seen  that  Eq.  (7)  is  verified. 

The  decorrelation  time  of  the  unsmoothed  temperature  difference 
fluctuations  is  fundamentally  on  the  order  of  the  (probe  spacing  ■ v)  ar.d 
will  normally  be  small  compared  to  that  corresponding  with  the  outer  scale. 

7.  S.  A.  Collins  and  G.  W.  Reinhardt,  "Investigation  of  Laser  Propagation 
Phenomena",  RADC-TR-7 1-248 , August  1971,  Electrrscience  Laboratory, 
Ohio  Ste.te  University.  (AD734  547) 
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Squared  temperature -difference  fluctuations  for  various  smoothing  times  (T). 

a.  t = 0.  4 sec. 

b.  t = 3. 2 sec. 

c.  T = 12. 8 sec. 
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ation  function  of  smoothed  fluctuations  as  in  Figure  6,  for  various 
times. 


* I 


Figure  8.  Computed  integral  scales  IT  of  the  smoothed  processes  corresponding 
to  Figures  6 and  7,  vs.  smoothing  time  (T  ). 
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Therefore,  in  practical  measurements  of  Cn  , the  integral  scale  will  be 
t/2.  The  same  relationship  between  the  averaging  time  and  integral  scale 
of  a short-term  average  will  obviously  also  apply  to  propagation  quantities 
such  as  the  log  amplitude  variance  of  scintillation,  which  tends  to  have 
an  integral  scale  which  is  an  order  of  magnitude  smaller  than  that  of  the 
thermal  fluctuations.  For  the  data  of  Fig.  8,  the  unsmoothed  microthermal 
scale  was  s 0.6  sec.,  while  the  optical  and  (10  micron)  infrared  scales 

were  40  and  60  msec  respectively. 

C.  Large-Scale  Intermittent 

The  practical  defin  tier,  of  an  ''intermittent"  or  "sporadic"  random 
process  is  that  such  a process  has  an  envelope  or  strength  which  varies 
much  more  slowly  than  the  process  itself,  but  changes  significantly  over 
observation  times  of  interest.  Since  the  concept  of  varying  "strength" 
or  "envelope"  implies  short-term  averaging  of  the  second  moment,  we  may 
relate  this  concept  to  the  relative  time  behavior  of  the  (short-term- 

averaged)  first  and  second  moments. 

The  question  of  long-term  stationarity  of  a process  is  often  unan- 
swerable and  moot;  for  mathematical  convenience,  we  assume  ultimate 
stationarity  and  model  the  process  in  terms  of  what  we  learn  during  the 
limited  period  of  observation.  It  is  apparent  that  intermittency  or 
sporadicity  as  defined  here  is  not  basically  related  to  questions  of  sta- 
tionarity. 

Atmospheric  turbulence,  as  manifested  in  the  temperature  difference 
fluctuations  between  two  probes,  often  exhibits  intermittency,  e.g.,  in 
the  presence  of  plumes,  or  in  poorly-developed  turbulence.  An  example 
is  shown  in  Figure  9a.  We  are  interested  in  applying  all  of  the  statistical 
considerations  of  this  discussion  to  the  intermittent  case  as  well  as  to 

the  simpler  condition  of  large-scale  uniformity. 

The  intermittent  process  has  what  at  first  are  discerned  to  be 
peculiar  properties.  The  slow  envelope  fluctuations  do  not  represent 
additive  low  frequencies  in  the  process,  but  rather,  are  indicati 
narrow  sidebands  around  high-frequency  spectral  components.  (Additive 
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low  frequencies  may  exist  in  single-probe  microthermal  measurements,  and 
may  relate  to  stationarity  problems,  but  they  do  not  relate  to  the  scintil- 
lation phenomena  of  interest  here.)  Since  low  frequencies  are  not  present, 
the  autocorrelation  function  of  the  process  does  not  have  a tail  and  does 
not  manifest  the  intermittency.  Hence,  as  will  be  discussed  in  a later 
section,  the  averaging -time  considerations  of  Ref.  7 must  be  used  with 
care  in  such  a case. 

In  order  :o  exhibit  the  low  frequencies  in  the  envelope,  we  must  per- 
form a nonlinear  operation,  such  as  taking  the  absolute  value  (Fig.  9b)  or 
squaring.  This  effectively  convolves  the  narrow  sidebands  around  each 
high-frequency  component,  resulting  in  additive  low  frequencies.  The 
autocorrelation  function  then  exhibits  the  effects  of  the  intermittency. 

As  discussed  in  the  following  section,  an  intermittent  process  (z)  of  this 
type  may  be  modelled  as  the  product  of  (xy)  of  fast  (x)  and  slow  (y) 
processes,  where  the  latter  (envelope)  process  is  always  zero  or  positive. 
This  model  may  in  turn  be  related  to  the  integral  scales  of  interest. 

In  the  simplest  model  for  intermittent  turbulence,  the  slow  envelope 

process  of  the  preceding  paragraph  has  only  two  values,  0 and  1.  Although 

0 

this  has  on  occasion  given  a reasonable  fit  to  experimental  data,  we  find 
that  in  general  it  is  not  a good  approximation.  To  explore  this  question 
further,  we  consider  a typical  example.  In  Fig.  10a  we  again  show  the 
fluctuations  of  Fig.  9,  and  the  absolute  value  smoothed  by  means  of  digital 
RC  filtering  with  progressively  longer  averaging  times.  The  same  is 
shown  in  Fig.  10b,  using  digital  boxcar  integration  for  smoothing.  In 
order  to  avoid  smearing  the  character  of  the  intermittent  envelope,  the 
maximum  smoothing  time  shown  is  0.5  sec. 

The  probability  distributions  for  the  basic  and  smoothed  processes 
of  Fig.  10  are  shown  in  Fig.  11.  The  distribution  function  for  the 
envelope,  with  0.5  sec  smoothing,  is  seen  to  have  a slight  tendency  for  a 
double-humped  appearance,  but  is  poorly  approximated  by  a two-level 
function. 


8.  R.  S.  Lawrence,  G.  R.  Ochs,  S.  F.  Clifford,  J.  Opt.  Soc.  Am.  60,  826 
(1970)  . 
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= 5 msec 


= 50  msec 


= 500  msec 


Figure  10a. 


iLa ^ 


Fluctuations  of  Fig.  9b.  for  increasing  smoothing  times  (t). 

a.  Smoothed  with  digital  RC  filter. 

b.  Smoothed  wiai  digital  boxcar  integrator 


Figure  10b. 


If  we  again  consider  the  multiplicative  model,  we  postulate  that  the 
fast  process  <x)  has  the  high-flatness-factor  statistics  of  the  fundamental, 
unsmoothed  microthermal  process  of  Sec.  II-A.  The  envelope  distrlhution  (y) 
will  be  maximum  at  small  values  and  will  have  a distribution  like  that 
<0.5  sec,  in  Fig.  11.  The  probability  distribution  (P,)  of  the  composite 

process  (z  = xy)  may  be  written  as 


P dz 
z 


= f p (z/y)  p (y)  1/y  dy  dz 
J X y 


(8) 


We  therefore^pect  the  composite  Intermittent  process  to  have  a probability 
distribution  with  a fairly  high  flatness  factor  but,  owing  to  the  envelope 
distribution,  it  will  have  more  smearing  or  less  of  a log  normal  li  e 

character  than  in  the  non-intermittent  case. 

Although  we  do  not  yet  have  a sufficiently  long  data  record  of  in  e - 
mittent  fluctuations  to  generate  a beta  diagram  vs.  smoothing  time  for 
squared  fluctuations,  such  as  that  of  Fig.  4 for  the  non-intermittent  case, 
we  can  speculate  on  the  results  of  such  a plot.  The  unsmoothed  fluctuations 
will  have  a smaller  flatness  factor  than  in  the  non-intermittent  case. 

For  intermediate  smoothing  times  (such  as  0.5  sec  in  Fig.  10),  which  exhibit 
the  intermittent  envelope,  the  same  will  be  true.  However,  as  discussed  in 
the  following  section,  the  low  frequencies  relating  to  intermittency  w 
cause  the  integral  scale  of  these  squared  fluctuations  to  be  larger  than 
in  the  non-intermittent  case,  and  hence  the  approach  to  a normal  distr  u 
tion  with  increased  smoothing  times  will  be  more  gradual.  These  compar  sons 
of  course  assume  an  identical  wind-speed  normalizing  factor. 

We  note  that  the  meteorological  and  fluid  mechanical  literature  con- 
tains many  discussions  of  microscale  intermittency  and,  in  some  cases  its 
relationship  to  large-scale  intermittency.  Many  of  these  have  been  refer- 
enced in  earlier  reports  on  this  program.  The  field  is  clearly  in  a state 
of  flux,  and  there  is  much  room  for  more  fundamental  work.  Our  approach  to 
macroscale  intermittency  is  a pragmatic  one,  in  that  it  avoids  esoter  c 


9.  G.  J.  Hahn  and  S.  S.  Shapiro,  Statistical  Models  in  Engineering, 
Wiley  & Sons,  1967. 
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questions  of  generalized  stationary , the  inadequacies  of  the  energy- 
cascade  model,  etc.10  A basic  complication  is  that  the  behavior  of  inter- 
mittent microthermal  envelopes  cannot  be  expected  to  be  universal  but 
rather,  will  be  a complicated  function  of  meteorological  and  local  geo- 
graphical factors.  Similarly,  the  localized  regions  of  stronger  or 
weaker  turbulence  may  not  be  expected  to  be  isotropic. 

Contrary  to  some  statements  in  the  literature,  it  now  seems  clear 
that  intermittency  plays  no  fundamental  role  in  such  propagation  phenomena 
as  the  saturation  of  scintillations  or  the  generation  of  log  normal  scin- 
tillation statistics.  We  may  therefore  view  the  role  of  intermittency 
as  simply  a mechanism  which  increases  the  fluctuations  in  short-term- 
averaged  propagation  measurements  or  channel  performance,  and  this  can  be 
quantitatively  described  in  terms  of  relevant  measures  such  as  spatial 
envelope  correlations  and  integral  scales.  These  factors  will  be  discussed 

further  in  following  sections. 

D.  Further  Discussion  of  Integral  Scales 

In  a previous  section  we  have  shown  that  the  integral  scale  of  any 
smoothed  process  approaches  one-half  the  averaging  time  as  the  latter 
increases.  We  have  also  pointed  out  that  a slowly  varying  envelope  or 
strength  of  a process  does  not  have  an  effect  on  the  autocorrelation,  low- 
frequency  content,  or  integral  scale  of  that  process.  However,  once  we 
take  the  absolute  value  (rectify)  or  square  the  process,  the  low  envelope 
frequencies  will  be  directly  manifested.  An  example  of  the  enhanced  low 
frequencies  in  the  absolute  value  of  intermittent  microthermal  fluctuations 

is  shown  in  Fig.  12. 

As  a prelude  to  the  next  section,  in  which  averaging  time  effects  in 
the  presence  of  intermittency  are  discussed,  we  show  here  analytically 
how  the  integral  scale  of  an  intermittent  process  is  affected  by  squaring. 
We  again  model  the  process  as  the  product  of  a fast  process  (x)  and  slow 
envelope  process  (y) . Since  these  processes  are  independent,  and  x has 
zero  mean,  the  composite  process  z = xy  also  has  zero  mean. 


10.  Benoit  Mandelbrot,  "Intermittent  Turbulence  in  Self  Similar  Cascades 
J.  Fluid  Mechanics,  62 .,  Part  II,  331-358  (1974). 


-24- 


trum  of  temperatur 


To  find  the  integral  scale  ot  z,  we  write 

f (x (t+t ' ) x (t ) y ( t+t 1 ) y (t)^  dt' 

" (x2)  (y2) 

(y2)f  (x(t+t’)x(t))  dt' 


(9) 


This 


does  not  manifest  the  intermit tency. 


as  expected.  , uave 

L the  integral  scale  of  the  fluctt  lone  of  . . 

2 2 

f k2(t+t')x2(t)y2(t+t’)y2(t))  - (*2)  (y  ) 

FToTi 


dt' 


(10a) 


I 2 
z 


Noting  that  a and  y are  independent,  we  have 

2 2 

f |(y2(t«,)y2(t)>  (*2(t+f)*2(t>)  - (*2)  (>2) 

[^v<7T^  <*2> 


>dt' 


(10b) 


I 2 

X 


We  expand  this  to  write 


I 2 
z 


f\  (y2(t+t')y2(t))  (x2(t+f)x2(t)-  (x2) 
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V.  -f  ■ ' 


1 


+ / | (x2)  (y2(t+t’)y2(t)-  ( y 2)  ) |dt 1 


.</>(/)  - <*2>V)2J 


(ID 


= I 2 + I 2 

Z1  z2 


Let  us  now  denote  the  integral  scales  for  the  fluctuations  of  x and 

2 

y about  their  means  as  1^2  and  I 2 respectively.  1^2  will  be  on  the 

order  of  I , while  I 2 will  be  much  greater  since  y is  slowly  varying, 
x y 

Consider  the  first  term  of  Eq.  (11).  The  first  or  y-expectation  will  be 
^y^  for  all  t'  such  that  the  second  or  x-expectation  is  nonzero.  We 
thus  have 


1 2 . (y4)  ’ 1xz(<«4)  - (*2 ) ) 

21  <*4>  <y‘>  - <«2)2  (y2>2 


I 2 
x 


1- 


<# 


- 1 


<v2)2 


Ix2  (1  - a) , 


(12a) 


where  1®  the  flatness  factor  defined  previously.*  For  the  second 

term  of  Eq.  (11) , we  have 


*The  term  ^y^  / ^y2^  is  not  identical  to  • since  ^y^  4 0.  It 


may 


be  taken  as  a measure  of  intermittency. 
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(12b) 


To  examine  the  behavior  of  these  expressions,  we  note  that  >:>  1* 
while  ^ 'y = 1 for  no  intermittency  and  is  greater  in  £he  pres- 
ence of  intermit tency.  Hence  for  no  intermittency,  a = 0 and  I 2 = I 2 = 

/ A _ z *- 1 

Ix2  as  expected.  As  \^L  grows,  I*  becomes  somewhat  less  than  I 2, 

<>2> 

while  I 2 grows  to  a maximum  of  I 2/6-  • Hence  the  integral  scale  or 
z y 2 

2 x 

decorrelation  time  of  an  unsmoothed  intermittent  process  when  squared  will 
exhibit  the  long  integral  scale  of  the  envelope,  diminished  by  the  high 
flatness  factor  of  the  basic  process.  This  will  be  pertinent  in  the  con- 
siderations of  later  sections. 

Finally,  we  point  out  that  the  behavior  of  the  integral  scale  when  a 
random  process  is  squared  can  serve  as  a quantitative  indication  of  the 
degree  of  "intermittency"  of  the  process;  Iz2  depends  on  both  the  envelope 

I 

z 

distribution  (a)  and  on  the  envelope  scal<  (1^2). 

E.  Averaging  Time  vs.  Data  Spread  for  Intermittent  Turbulence 
Considerations  of  averaging  time  effects  on  measurements  of  random 
processes  have  been  reviewed  in  Ref.  7.  At  first  appearance,  there  has 
been  a paradox  in  the  application  of  these  results  to  intermittent  processes, 
in  that  data-spread  is  expressed  in  terms  of  the  integral  scale  (or  power 
spectrum)  of  the  process,  while  the  envelope-intermittency  does  not  show 


/ 
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up  in  these  quantities.  Yet,  it  is  physically  clear  that  intennittency 
can  strongly  affect  the  spread  in  measurements  of  short-term-average 
quantities,  even  though  the  process  is  truly  stationary  in  the  long  term. 
The  purpose  of  this  section  is  to  resolve  this  seeming  inconsistency,  and 
to  remark  on  the  interpretation  of  the  results  for  this  case. 

Let  us  consider  the  stationary,  intermittent  process  z(t)^_We  denote 
any  sample  measurement  of  zn(t)  over  an  averaging  time  x by  z"  . The 
integral  scale  of  the  fluctuations  of  zn  is  again  denoted  by  I n]  and  we 

assume  that  any  such  integral  scale  is  short  relative  to  averaging  times 
of  interest. 

If  we  now  define  f(t)  in  Ref.  7 as  z(t),  and  recall  that  (z).  0, 
a simple  statement  of  the  results  of  Eq.  (8)  of  that  reference  is  that 


<4> 

<°2>  ' T 

This  is  not  a very  interesting  result,  because  the  simple  deviation  of 
sample  averages  of  z(t)  from  its  true  mean  of  zero  is  not  interesting. 
The  dependence  on  ^ which  does  not  manifest  the  intermittency,  is 
clearly  appropriate  in  the  case  of  this  measurement. 


(13 


that 


We  now  redefine  f(t)  of  Ref  7 a<?  ,2/v\  „ 

vw  ui.  Ret.  / as  z (t).  We  then  have  the  result 


<( 


(a4)  * 


<>2)2 


2 12 
z 


and' 


Var 


<’2> 


2 1?. 

z 


Var  : 

<‘2> 


(14) 
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This  is  the  result  of  interest,  and  as  discussed  in  the  previous  section, 

I 2 does  manifest  the  intermittent  envelope.  To  apply  this  result,  we 
z 

must  know  the  long-term  or  approximate-ensemble  values  involved.  For  an 
intermittent  process  with  slow  envelope  variations  (large  I 2),  this  may 
require  a very  long  time.  Also,  the  ensemble  quantity  / 2\  does  not 


<■/> 


apply  to  a limited  set  of  measurements,  and  in  the  presence  of  slow  inter- 
mit tency,  may  likewise  require  many  samples  to  determine  approximately. 
However,  there  is  no  formal  inconsistency. 

We  point  out  that  the  self-consistent  averaging -time  argument  in 
Ref.  7 implicitly  assumes  that  no  long-term  variations  in  the  envelope 
of  z(t)  occur.  I.e.,  within  a time  span  shorter  than  the  long-term  fluc- 
tuations in  the  envelope  fluctuations,  whose  existence  might  not  even  be 
known  to  the  experimenter,  the  ensemble  quantities  are  an  abstraction; 
there  is  no  way  to  know  their  true  values. 

We  also  point  out  that  the  condition  that  t is  much  greater  than 
the  relevant  integral  scale  is  much  more  severe  for  Eq.  (14)  than  for  (13). 
However,  shorter  averaging  times  (e.g.  corresponding  to  the  outer  scale  of 
turbulence)  may  retain  more  information  and  therefore  be  of  interest. 

Hence,  in  the  presence  of  intermittency,  the  averaging  time  may  not  con- 
tain many  decorrelation  times  for  the  envelope.  An  appropriate  generaliza- 
tion of  the  above  results,  not  requiring  t » I,  may  be  written  from 
Ref.  7: 


(i\  -<z2))  ) 2 

9 x 

“ T 

f (1-  t/x)  p 2 (t)  dt 

44>  42> 

L-l  z J 

(14a) 


2 ? 

where  p 2 is  the  autocorrelation  of  z . In  the  limit  t <<  14,  this 
z z 

equals  Pz2  (0),  which  will  be  much  larger  than  the  1^2  « T result  of 
Eq.  (14). 
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F.  Fourier-Stieljes  Representation  for  Intermittent  Turbulence 
As  mentioned  above,  our  operational  description  of  intermittency 
circumvents  the  need  for  the  complications  of  generalized  stochastic 
theory.6,  10,  11*  12*  Although  it  was  not  originally  obvious,  it  may  be 
easily  shown  that  intermittency,  as  modelled  here,  does  not  introduce 
difficulties  in  the  conventional  spectral  representation  for  a random 
process.  This  may  be  shown  as  follows. 

The  conventional  Fourier-Stieljes  spectral  representation  for  a 

stationary  random  process  f(t)  is  given  by 

TO 

f(t)  = j e iWtdF(u>) 

— 00 

where  dF(w)  is  a random  function  of  the  increment  in  w.  More  explicitly, 

or  Q 

e^UtdF((i)) 

and 

n N 

( e^'dFM  = lim  V e^k'1  [ F(^H(Vl)l  (15) 

in  | |a»kl  1-0 

where 

-fi  = 0JQ<a)1<a)2 

“A'^-l 

| | Aiok|  | = max  l“k-U)k_1l 

Since  the  autocorrelation  of  f(t)  cannot  depend  on  t,  the  random  variables 
F(uj)  and  F((d._1)  must  be.  uncorrelated  for  This  rollows  from 

11.  A.  M.  Yaglom,  An  Introduction  to  the  Theory  of  Stationary  Random 
Functions,  Prentice-Hall,  Inc. , 1962. 

12.  J.L.Doob,  Stochastic  Processes,  Wiley  & Sons,  1953. 
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writing  the  autocorrelation  as 


CO  00 

(t)f*(t+t')^  ~/f  e1U)tdF(o))  y*e"lw,(t+t,)  dF(w') 

'—■00  —00 


00 


e i^t'git^  w ' ) ^dF (w)dF (w' ))  . (16) 

—00 


Hence,  we  have  the  familiar  requirement  that  the  spectrum  of  a (wide- 
sense)  stationary  process  have  "orthogonal  increments": 

^dF(o))dF((jj' = 6(u)-u)')^|dF(o))  | ^ . (17) 

Now  consider  the  multiplicative  model  of  an  intermittent  process, 
z = xy.  Since  x(t)  and  y(t)  are  independent  stationary  processes,  the 
autocorrelation  of  z(t)  may  be  written  as 

^z(t)z(t+t')^  = ^x(t)x(t+t')^  ^y(t)y(t+t')^ 

= Px(t’)Py(t’) 

= p (t ' ) . (18) 

z 

Hence  z(t)  is  trivially  shown  to  be  stationary,  with  a spectrum  which 
also  has  orthogonal  increments. 

The  key  to  this  result  is,  of  course,  the  independence  of  the  fast 
fluctuations  and  the  modulating  envelope.  If  they  are  stationary  but  not 
independent,  the  probability  distribution  of  the  product  (Eq.  8)  involves 
a conditional  probability  which  in  turn  involves  the  joint  distribution  of 
x and  y.  It  is  known  that  the  latter  may  not  be  stationary  in  such  a case, 
so  that  z = xy  may  not  be  stationary. 


13.  A.  Papoulis,  Probability  Random  Variables  and  Stochastic  Processes, 
McGraw-Hill,  1965. 
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Most  of  our  remarks  here  involve  double-probe  (or  high-passed) 
microthermal  data  appropriate  to  scintillations,  in  which  independence 
of  the  envelope  and  basic  fluctuations  is  probably  a valid  assumption. 
However,  in  the  case  of  single-probe  data,  which  relates  to  optical/ 
infrared  phase  fluctuations  and  which  has  additive  low  frequencies, 
this  independence  is  clearly  violated,  e.g.,  in  the  case  of  plumes  of 
warmer,  more  turbulent  air  relative  to  the  surroundings. 

We  intend  to  subject  the  hypothesis  of  orthogonal  increments  to 
computer  testing,  for  single  and  double-probe  fluctuations  exhibiting 

interraittency. 

G.  Spatial  Correlations  of  the  Microthermal  Envelope 

As  will  be  seen  in  a later  section,  the  most  important  property  of 
the  spatial  microthermal  envelope  is  its  two-point  correlation  function 
down  the  propagation  path.  This  determines  the  number  of  independent 
"cells"  traversed  by  the  path,  and  consequently  the  fluctuation  in  the 
short-term  strength  of  scintillations.  Since  measurements  of  the  micro- 
thermal  field  per  se  are  normally  made  with  instrumentation  at  one  or 
more  points,  where  the  wind  transports  the  turbulence  structure  past 
each  probe,  we  are  also  interested  in  the  "frozen-in"  or  Taylor  hypothesis 
philosophy  as  it  applies  to  the  envelope.  Similarly,  since  the  wind 
direction  is  in  general  randomly  related  to  the  path  vector,  we  are 
interested  in  the  horizontal  isotropy  of  the  envelope. 

These  considerations  of  isotropy  and  taylor  hypothesis  as  applied 
to  the  microthermal  envelope  are  quite  different  from  those  for  the  basic 
fluctuations  themselves.  The  term  "envelope"  implies  some  degree  of  time 
or  spatial  averaging,  and  the  temporal-spatial  scales  involved  are  larger 
than  the  outer  scale  <LQ) . Hence  fundamental  similarity  arguments  are 
not  applicable,  and  universal  behavior  can  hardly  be  expected. 

Consider  the  microthermal  sensor  configuration  of  Fig.  13a,  where 
each  point  represents  a two-meter  high  differential  probe  pair,  for  which 
the  output  is  squared  and  short-term-averaged.  This  squared,  smoothed 
output  is  then  cross-correlated  with  that  from  each  other  point,  as  a 
function  of  time  lag.  When  little  or  no  smoothing  is  used,  and  the  point 
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separations  (AL)  in  each  direction  are  1.5  «,  the  normalized  cross  corre- 
lation for  adjacent  points  (12  or  23)  is  typically  0.3  at  zero  time  lag; 
it  falls  monotonically  for  points  23  perpendicular  to  the  wind,  while 
for  points  12  along  the  wind,  it  has  the  same  zero-lag  value  and  exhibits 
a secondary  peak  ,f  typically  0.1  at  a time  lag  corresponding  roughly  to 
the  "Taylor  interval"  (separation  divided  by  wind  speed).  However,  if 
the  point  separation  is  increased  to  2.5  m,  the  zero-lag  correlation  is 
very  small  (0.02)  and  independence  is  exhibited  at  all  time  lags.  Hence, 
for  the  squared  but  unsmoothed  case,  where  the  additive  low  envelope  fre 
quencies  are  dominated  by  the  high  frequencies,  a complete  breakdown  in 
correlation  is  seen  for  separations  on  the  order  of  the  outer  scale.  In- 
side this  scale,  very  good  isotropy  and  significant  frozen-in  structure 

are  observed. 

Now  consider  the  envelope  correlations  with  significant  averaging 
time.  Just  as  the  averaging  time  affects  and  can  determine  the  integral 
scale  for  an  autocorrelation  function,  it  may  also  be  expected  to  affect 
the  isotropy  and  frozen-in  character  of  the  spatial  envelope.  Consider 
the  probe-pair  measurement  points,  2 meters  high  and  separated  as  in 
Figure  13b  (eight  probes  total).  This  arrangement  provides  point  separa- 
tions of  2.5,  4,  6,  6.5,  10,  and  12.5  m.  Using  0.4  sec  averaging  with  a 
3.5  m/sec  wind  along  the  line  joining  the  points,  the  respective  mean- 
square  fluctuations  were  as  shown  in  Fig.  14.  The  Taylor  time  lag  and 
decaying  correlation  are  qualitatively  evident.  Note  that  the  averaging 
time  corresponds  spatially  to  a scale  (1.4  m)  on  the  order  of  Lq. 

The  time-lagged,  spatial  correlation  functions  are  shown  for  shorter 
separations  in  Fig.  15a  and  for  larger  separations  in  Fig.  15b.  The 
spatial  and  temporal  autocorrelation  functions  are  shown  in  Fig.  1* . 
showing  a breakdown  in  agi  ement  at  approximately  8 m.  Clearly,  smoothing 
increases  the  realm  of  frozen-in  envelope  behavior  substantially,  by 
emphasizing  the  low  frequencies  in  the  envelope.  The  spatial  vector- 
correlation  effects  (isotropy)  of  averaging  will  be  considered  further  in 
the  next  section.  The  details,  of  course,  will  depend  on  the  behavior  of 
the  (two-dimensional)  spatial  frequencies  below  that  determined  by  the 
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Figure  15a. 

Time -lagged  spatial  correlation  functions  C(t  ) for  fluctuations  of 
Fig.  14. 


a.  Smaller  measurement-point  separations 

b.  Larger  separations. 


smoothing  time  and  wind  speed;  in  the  presence  of  pronounced  intermit- 
tency,  still  larger  scales  will  predominate.  Further  two-dimensional 
experiments  are  being  conducted. 

H.  Microthermal  vs.  Propagation  Data  Spread 

In  this  section  ve  will  directly  address  the  problem  of  data  spread 
for  finite  averaging  times,  for  both  microthermal  measurements  of  turbulence 
strength  (C  2)  and  optical/infrared  measurements  of  scintillation  strength 
(a2  = log  amplitude  variance),  and  their  interrelationship.  We  begin  with 
purely  empirical  data,  and  then  we  present  an  analytical  treatment  along 
with  preliminary  data  of  a more  sophisticated  nature. 

When  reasonably  long  averaging  times  are  used  relative  to  the  integral 
scales  of  the  processes  involved,  one  would  hope  for  good  agreement  between 
theoretical  predictions  and  exper imetital  results.  In  Fig.  17,  we  show 
the  measured  log  amplitude  variance  of  scintillation  vs.  the  value  as 

predicted  by  first-order  theory  for  a uniform  path: 


0.124 


.7/6  .11/6 
k L 


(19) 


where  k is  the  optical/infrared  wavenumber,  L is  the  pathlength,  and 

C 2 is  determined  from  microthermal  measurements  at  a point.  The  averaging 

n 7 i n 

time  used  for  determining  the  abscissa  and  ordinate  of  each  point  was  7-10 

minutes,  and  each  two-wavelength  data  run  was  performed  on  a different  day. 

The  points  for  10.6  microns  are  (with  one  exception)  free  from  the  effects 

of  saturation  of  scintillations,14  and  very  good  agreement  with  theory  is 

observed.  These  data  were  not  specially  chosen  in  any  way,  and  the  only 

requirement  on  the  meteorological  conditions  was  that  the  sky  be  clear  or 

overcast  in  a nominally  uniform  manner.  The  points  at  4880X  are  of  course 

indicative  of  saturation  at  the  1.6  km  pathlength  used. 

In  Fig.  18  we  show  a plot  of  most  of  the  same  10.6  micron  data  of  Fig. 17, 

with  an  indication  of  the  spread  observed  in  short-time-averaged  (10  sec) 

2 2 2 

measurements  of  both  o and  C or  o_  . The  bars  indicate  one  standard 

E n i 


14.  R.  S.  Lawrence  and  J.  W.  Strohbehn,  Proc.  IEEE  58,  1523  (1970). 
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deviation  in  a large  number  of  independent  10  second  measurements,  crossing 
at  the  mean.  We  observe  a much  greater  spread  in  the  abscissa  than  the 
ordinate,  as  expected,  since  the  former  relates  tc  a microthermal  measure- 
ment at  a point,  while  the  latter  relates  to  a propagating  beam  which 
experiences  a weighted  spatial  average  over  the  randomly  varying  turbulence 
strength  along  the  path.  This  is  also  manifested  in  the  integral  scale  of 
the  squared  optiral/inf rared  amplitude  (the  intensity)  being  much  smaller 
than  that  for  the  squared  thermal  fluctuations.  Ordinate  spreads  in  the 

saturated  4880&  regime  are  comparable  to  the  unsaturated  10.6  micron  spreads. 

2 

If  we  consider  the  abscissa  or  10-second  C measurements  alone,  for 

n 

this  fixed  averaging  time  we  would  expect  the  spread  to  tend  to  be  inversely 
proportional  to  wind  speed  (v) , since  v determines  the  time  scale  of  all 
of  the  fluctuation1’.  This  tendency  is  supported  by  the  data  of  Fig.  19, 
although  the  degree  of  variability  or  intermit tency,  and  its  scale,  are 
important  added  factors. 

Let  us  now  consider  the  analytical  relationship  between  expected 

microthermal  and  optical  data  spread.  We  rewrite  Eq.  (19)  in  a form  that 

2 2 

explicitly  indicates  the  varying  nature  of  and  o , for  any  specific 
realization  of  the  path: 

L 

* j f (z)Cjj2  (z)  dz  (20) 

o 

where  f(z)  is  a path-weighting  function  dependent  on  optical  geometry, 

and  t is  the  finite  averaging  time  of  the  measurements.  The  viewpoint 

here  is  that  the  path  is  statistically  homogeneous,  and  the  variations 

indicated  in  Eq.  (20)  are  simply  due  to  the  short  term  nature  of  the 

averages  involved.  We  also  assume  stationarity  and  ergodicity.  For  a 

spherical  wave,  f(z)  ^ (z/L)5^(L-z)5^. 

In  order  to  write  this  relationship,  it  must  be  assumed^  that  the 
2 

turbulence  strength  varies  slowly  over  a scale  size  on  the  order  of 

15.  V.  I.  Tatarski,  Propagation  of  Waves  in  a Turbulent  Atmosphere,  Nauka, 
Moscow,  1967. 
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the  outer  scale  L , and  it  is  implicit  that  the  mean-square  microthermal 
° ~ 

fluctuations  (or  C Z)  are  z-averaged  over  at  least  a comparable  scale 
n 2 
(t>l  /v).  If  we  think  of  Eq.  (20)  as  expressing  the  time  function  (t), 

— o 

it  is  also  clear  that  this  function  will  not  vary  significantly  over  a 
time  scale  shorter  than  L /v,  since  the  integral  does  not  contain  information 
on  faster  fluctuations  of  the  scintillation  strength.  Of  course,  t may 
be  arbitrarily  longer  than  this  minimum  value,  and  as  it  becomes  very 
large,  Eq.  (20)  becomes 

L 

f(z)dz  . (20a) 


lim 


0?  - <°2>  - (Of 


However,  this  loses  the  information  of  interest  in  this  discussion. 

2 

Let  us  now  consider  a series  of  repeated  measurements  of  and 

C 2,  the  latter  at  any  point  on  the  statistically  homogeneous  path,  with 
finite  t.  The  ensemble  mean  of  a large  number  of  such  measurements  will 
yield  the  same  result  as  Eq.  (20a).  Of  interest  is  the  relationship  between 
the  normalized  variance  (data  spread)  in  the  values  obtained  for  these  two 
parameters.  We  therefore  write 


It  is  apparent  that  the  data  spread  will  be  determined  by  the  spatial  cor 
relation  function  of  C 2 vs.  z,  for  the  averaging-time  specified. 

We  may  also  rewrite  Eq.  (21)  in  the  form 
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Var  C 
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n 
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(22) 


A useful  (discrete)  approximation  is  to  consider  the  z-correlation  of 
to  be  either  total  or  zero,  depending  on  whether  Iz^-z^  t^ie  corre^at^on 
length,  so  that  the  overall  correlation  term  on  the  right  side  of  Eq.  (22) 
is  either  unity  or  zero  respectively.  If  we  then  replace  the  weighting 
function  integrals  with  corresponding  discrete  summations,  we  have 


(23) 


where  N is  (1/2)  the  number  of  correlation  lengths  of  C in  L.  For  a 

nT 

uniform  path-weighting  function,  this  takes  the  particularly  simple  form: 


,,  2 
Var  a 

T 


Var  C 2 
n 

i 


1 

N 


(24) 


For  a general  weighting  function,  we  may  define  the  ratio  of  summations  in 

Eq.  (23)  as  N (<N). 

eq  — 
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In  order  to  discuss  the  effects  of  the  averaging  time  on  the  optical 

data  spread,  and  the  relation  to  pertinent  integral  scales,  we  first  note 

2 2 

that  the  fundamental  random  variables  related  to  C and  a are  the 

0 n t 

L.  T 

squared  temperature  fluctuations  (AT)  and  the  squared  log  amplitude 

fluctuations  (log  A - (log  a)  ) respectively.  The  integral  scale  (1^2) 

of  the  former  may  he  >>L  , depending  on  the  degree  of  intermittency , and 

o 

the  integral  scale  of  the  latter  will  depend  on  the  wind  direction.  We 

now  distinguish  two  cases: 

1.  Wind  perpendicular  to  the  path 

In  this  case,  the  z-correlation  and  hence  N,  is  unaffected  by 
the  averaging  time  x.  The  integral  scale  of  the  squared 
optical  fluctuations  will  be  1^2,  and  for  x>I^T2,  both 
normalized  ensemble  variances  or  data  spreads  in  Eqs.  (22-24) 
will  be  'v  x ^ in  accordance  with  Eq.  (14). 

There  is  a subtlety  leading  to  an  apparent  contradiction 
here.  When  we  specify  x > 1^2,  we  are  specifying  that  the 
averaging  exceeds  intervals  over  which  large-scale  or  inter- 
mittent variations  in  the  turbulence  strength  may  occur.  As 
the  averaging  time  is  increased,  in  accordance  with  Fig.  5 
progressively  lower  frequencies,  along  the  transverse  vector, 
are  emphasized.  One  may  expect  these  frequencies  or  long 
spatial  wavelengths  to  possess  some  quasi-isotropic  properties, 
i.e.  to  have  enlarged  perpendicular-  or  z-correlations.  This 
would  imply  that  N is  proportional  to  x \ so  that  optical 
data  spread  is  independent  of  averaging  time.  However,  these 
residual,  low-frequency  components,  even  though  emphasized 
owing  to  the  smoothing  process  (x),  are  not  correlated  along 
a direction  (z)  perpendicular  to  the  smoothing. 

2.  Wind  along  the  path 

In  this  case,  if  L > x > 12 > then  in  accordance  with  the 
discussion  of  SecY  B,  the  z-correlation  is  equal  to  vx/2. 

This  assumes  that  the  turbulence  strength  or  envelope  satisfies 
the  "frozen-in"  viewpoint.  From  Eq.  (14),  N = L/vx,  and 
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2 

Var  o 2/{rj2)  wil1  be  independent  of  averaging  time.  This  manifests  the 
fact  that  the  integral  scale  of  the  squared  optical  fluctuations  is  on 
the  order  of  L/v. 

When  x>L/v,  we  have  N=l,  and  the  optical  spread  will  be  . 

These  new  understandings  of  expected  optical  data  spread  are  appro- 
priate for  experimental  confirmation.  The  z-correlation  length  should  be 
directly  measured,  and  the  microthermal  and  optical  data  spreads  determined 
vs  t for  different  wind  directions.  Preliminary  results  will  be  dis- 
cussed below.  Unfortunately,  the  details  of  z-correlations  of  Cn  (z)  will 
depend  upon  wind  conditions,  terrain,  and  turbulence  and  hence  are  not 
universally  describable. 

It  is  interesting  to  point  out  that,  since  the  path-weighting  function 
[ f (z)  or  f.]  involves  wavelength  only  as  a multiplicative  constant,  the 
fractional  spread  in  log  amplitude  variance  measurements  are  predicted  to  be 
independent  of  wavelength.  Also,  note  that  the  representation  of  ^ as 

V*  f c ^ where  the  C 2 are  independent,  leads  to  the  calculation  of 

i-i  ^ °i  * ni  2 

the  moments  of  a2  distribution  from  those  of  the  C r distribution;  this 

T . T 
follows  from  the  characteristic  functions  <fi: 

N 

<f>  2 (s)  = TT  Vr  C (s) 

o i.  n- 

T i = 1 11 

where  *(s)  are  the  respective  characteristic  functions.  We  expect  that 
o'  tends  to  a gaussian  distribution  in  many  cases.  Finally,  the  develop- 
ment leading  to  Eqs.  (22-24)  is  also  applicable  to  the  case  of  saturated 
scintillations,  once  the  appropriate  weighting  function  (f)  is  known. 

It  is  apparent  that,  for  a general  wind  direction,  the  maximum  value 
of  N is  L/vzx,  and  when  vz  is  small  and/or  the  turbulence  strength  correla- 
tion regionZis  large,  N will  be  arbitrarily  smaller  than  this  value.  In 
Table  I,  we  show  preliminary  data  of  N^L/v^,  and  N approximated  from 
Eq.  (24).  The  wavelength  was  10.6  y,  the  pathlength  was  1600  meters,  and 

16.  H.  T.  Yura, "Ir radiance  Fluctuations  of  a Spherical  Wave  Propagating 
under  Saturation  Conditions",  to  be  published  in  JOSA. 
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t was  10  sec. 


TABLE  I 


Comparison  of  N determined  from  relative  optical/microthermal  data  spread 
in  repeated  measurements,  with  Nmax  determined  from  path-vector  component  of 
wind  speed. 


C 2 (m 
n 

-2/3) 

V (m/ 

Total v 

1. 3 x 

io'14 

5.8 

1. 8 x 

io-13 

0.8 

8.3  x 

10-14 

2.5 

X 

o 

00 

io-13 

1.8 

1 . 0 x 

io-13 

8.0 

2.9  x 

io-13 

6.0 

1.9  x 

io-12 

3.5 

3.9  x 

io-13 

1.8 

1. 5 x 

io-14 

6.0 

V (m/sec) 

N 

N 

z 

max 

5.4 

29.6 

40.0 

0.8 

20.0 

29.6 

2.24 

71.4 

22.8 

1.8 

88.9 

61.5 

7.8 

20.5 

24.2 

5.0 

32.0 

12.1 

3.2 

50.0 

33.0 

1.8 

88.9 

80.0 

5.4 

29.6 

19.0 

From  this  table  it  can  be  seen  that  the  present  discussion  is  well  supported 
qualitatively;  N exceeded  only  twice.  The  z-correlation  length  (=L/2N) 

varied  from  10  to  66  meters,  and  the  reduction  in  optical  data  spread  vs. 
microthermal  spread  ranged  correspondingly  from  80  to  12. 

It  is  obviously  desirable  to  have  direct  measurements  of  the  z-correla- 
tion  when  the  optical  and  microthermal  spreads  are  compared.  As  of  this 
writing,  we  have  one  reliable  data  sample,  as  given  in  Table  II.  The  data 
relate  to  the  signals  and  correlations  shown  in  Figs.  14-16. 


TABLE  II. 

Sample  comparison  of  theoretical  vs.  experimental  data  spread. 

Weather  conditions:  Clear,  T=17°C,  wind  along  path  at  3.5  m/sec 

C 2 = 1.5  x 10'14m"2/3,  X = 4880& 
n 

Four  probe  pairs  spaced  at  z=0,  2.5,  6.5,  12.5m 
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1024  data  points  sampled  at  0.4  sec  intervals  (6  min  50  sec  total) 

Relative  standard  deviation  in  measured  from  4 respective  probes=0.025 
Averaging  time:  0.4  sec 


, „ 2 Var  C 2 

Mpasured  Var  0 n 


(from  Eq.  22)  = 0.0115 

In  order  to  test  the  hypothesis  that  normalized  data  spreads  are 

independent  of  wavelength,  experiments  were  carried  out  over  simultaneous, 

coincident  paths  of  1600  m length,  at  4880X  and  1016  p.  Unfortunately, 

the  shorter -wavelength  scintillations  were  saturated  in  most  conditions. 

The  results  are  shown  in  Fig.  20. 

I.  Data  Spread  and  Intermittency  in  Imaging  Applications 

Titterton  ^ has  raised  the  question  of  "beam  breathing",  related  to 

18 

variations  in  the  coherence  parameter  p^.  In  our  phenomenological 
theory  of  beam-wave  effects,^  we  have  not  included  "breathing  and  see  no 
noticeable  deficiency.  However,  although  such  a variation  is  in  principle 
included  in  any  long-term  or  ensemble  formulation  of  propagation  problems, 
it  would  be  of  interest  to  consider  the  short-term  behavior  of  Pq  and  the 

17.  Paul  J.  Titterton,  "Beam  Spread  Statist:)’  < from  Reciprocity", 

Topical  Meeting  on  Optical  Propaga  ■>.'  .ough  Turbulence,  Boulder, 
Colorado,  July  1974. 

18.  J.R.Kerr , J.R.Dunphy,  "Propagation  of  Multiwavelength  Laser  Radiation 
Through  Atmospheric  Turbulence",  RADC-TR-74-183,  May  1974.  (AD783  277) 
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short-term  statistics  of  e.g.  image  resolution.  (By  "short-term"  here 

we  are  not  referring  to  quasi-irstantaneous  image  determinations,  in  which 

18 

image-dancing  is  resolved.) 

The  dynamic  behavior  of  pq,  including  its  probability  distribution, 

is  not  described  by  the  results  of  time-separated  measurements;  an  example 

of  the  latter  is  in  Ref.  19,  in  which  p is  seen  to  be  log  normally  dis- 

o 

tributed  on  a dav-to-day  sample  basis.  To  get  at  this  dynamic  behavior, 

. * • 18 
we  note  that  for  a point  source: 

L 

p "5/3  = 1.45k2  f C 2(z)(z/L)5/3dz  . (26) 

o In 

o 

We  may  therefore  use  the  methods  of  the  preceding  section,  but  the  results 

will  apply  to  the  (-5/3)  power  of  the  desired  variable.  As  a further  degree 

18 

of  sophistication,  one  may  utilize  the  so-called  "short-term  MTF",  which 
follows  the  centroid  of  a dancing  image  (i.e.  removes  the  linear  component 
of  phase  distortion).  This  topic  comprises  a possibility  for  future  work. 

HI.  RECENT  RESULTS  ON  LASER  ILLUMINATION  THROUGH  TURBULENCE 

18 

In  the  preceding  report,  a comprehensive  discussion  was  given  of  our 

current  understanding  of  turbulence  effects  on  finite  laser  beam  illumination 

systems,  including  wander-cancelling  systems.  Also,  the  first  set  of  results 

of  an  ongoing  experimental  effort  were  reported.  In  this  section,  two 

recent,  further  analytical  results  are  presented. 

A.  Effective  Aperture  for  a Truncated  Gaussian  Transmitter 

As  discussed  in  Ref.  18,  the  results  for  the  mean  on-axis  target  irra- 

diance  for  a focused  beam  propagating  through  turbulence,  with  and  without 

wander  cancellation,  can  be  expressed  universally  in  terms  of  the  ratio 

D/p  as  an  independent  variable.  The  parameter  D is  some  measure  of  the 
o 

transmitter  diameter,  and  pq  is  the  point-source  coherence  length  over 
the  reciprocal  path  (i.e.,  from  target  to  transmitter).  We  explicitly 
considered  both  a uniform  circular  aperture,  and  a nontruncated  gaussian 
aperture,  and  found  that  the  results  are  practically  the  same  in  both 
cases  once  the  proper  comparative  abscissas  (effective  diameters)  have 

19.  D.L. Fried  and  n.F,.Mevers,  "Evaluation  of  r for  Propagation  Down  Through 
the  Atmosphere",  preprint,  Nov.  1973. 
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been  established.  It  is  then  reasonable  to  expect  that  these  results  will 
also  apply  to  the  practical  case  of  a truncated  gaussian  transmitter;  the 
purpose  of  this  section  is  to  derive  the  appropriate  measure  of  diameter 
of  the  beam  for  the  general  case  of  arbitrary  truncation. 

To  do  this,  we  write 


1(a) 


dx  x f (x)e 


-(ax) 


5/3 


(27) 


1 R 

which  is  obtained  from  the  Huygens-Fresnel  formulation..  Our  notation, 
which  is  slightly  changed  from  Ref.  18,  is  as  follows: 

Dg  = 2a  = 1/e  diameter  of  gaussian  intensity  distribution  at  transmitter 
D = physical  diameter  of  truncation  aperture 

x = p/D  = normalized  difference-coordinate  in  transmitter  plane 

O 

a = D /p 
g o 

The  function  f(x)  is  from  the  overlap  integral  for  the  particular  distribution 
of  illumination  over  the  aperture.  Although  the  above  expression  could  be 
worked  out  rigorously  for  each  value  of  D/D  , the  near  coincidence  of  non- 

s 

truncated-gaussian  and  uniform-illumination  cases  suggests  that  we  need 

merely  define  an  appropriate  "effective  diameter"  for  the  general  case. 

To  do  this,  we  note  from  Ref.  18  that  the  effective  diameter  (D  ..) 

ef  f 

for  both  the  uniform  case  and  the  noniruncated  gaussian  case  is  specified 
by  the  ratio  of  the  limits: 


lim 

or*” 

lim 
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1(a) 
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where  Dgff=D  and  2D  for  the  two  cases  respectively.  For  the  general 

O 

case,  we  define 
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Then  from  the  truncated-beam-wave  diffraction  theory  in  a vacuum,  we 
write 


I 
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U k 
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(l-e_n/2)' 


(29) 


where  U is  the  transmitter  amplitude  on  beam  axis.  We  note  that,  from 
simple  integration,  the  transmitter  power  P is  given  by 


P * 
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For  the  nontruncated  case. 


we  can  write 
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We  know  that  I /P  is  controlled  by  the  turbulence,  and  thus  for  the  trun- 
cated case  this  ratio  must  be  independent  of  D,  D . It  is  thus  apparent 

8 

dimensionally  that  for  the  general  case. 
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(32) 
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As  a check,  we  note  that  this  approaches  the  correct  limit  when  tt*0. 


20.  G.  Olaofe,  J.  Opt.  Soc.  Am.  60,  1654  (1970). 
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We  now  take  the  Ratio  of  Eqs.  (32)  and  (29),  and  combine  with  (28), 
to  write 


eff 


2D 


1-e 


-n/2 


(33) 


The  ratio  Deff/2Dg  is  plotted  in  Fig.  21.  These  results  will  be  used  in 
the  future  to  unify  all  presentations  of  data.  They  may  also  be  used  as 
an  approximation  to  the  effective  aperture  for  a truncated  gaussian  beam, 
when  second  moment  or  irradiance  fading  data  are  presented. 

Finally,  we  note  that  the  phenomenological  theory  given  in  Ref.  18 
can  be  written  in  terms  of  rQ,Deff  instead  of  PQ,Dg,  with  a resulting 
simplification  in  the  numerical  constants.  This  will  be  done  in  a forth- 
coming publication  based  on  that  report. 

B-  Lfoit  of  Irradiance  Variance  at  Large  Values  of  D/o 

In  the  previous  report,18  it  was  pointed  out  that  both  the  phenomeno- 
logical model  and  experimental  data  suggest  that  the  normalized  irradiance 
variance  for  a transmitter  much  larger  than  the  coherence  scale  must  ap- 
proach unity.  Through  reciprocity,  this  would  also  be  true  of  a large 
heterodyne  reciever.  Although  this  seems  to  be  clearly  necessary  physically, 
it  disagrees  with  earlier  approximate  analyses  of  heterodyne  systems.21 
The  point  of  the  present  section  is  to  derive  this  limit  directly  from  the 
Huygens-Fresnel  expression  for  the  second  moment  of  irradiance  (fourth 
moment  of  amplitude) . 

The  required  Huygens-Fresnel  expression  is  given  in  Ref.  22.  For  a 

nontruncated  gaussian  beam,  this  may  be  readily  reduced  to  the  eight -fold 
integral: 


21. 


22. 


' ^r»ed>  Effects  of  Atmospheric  Turbulence  on  Static  and  Tracking 
Optical  Heterodyne  Receivers/Average  Gain  and  Antenna  Gain  Variation", 
echnicai  Report  No.  TR-027,  Optical  Science  Consultants,  August  1971. 
H.  T.  Yura,  Applied  Optics  11,  1399-1406,  June  1972. 


-55- 


(.’HA)  JJ(h 


d->-  d->  d->-  exp  < - 

P2  P3  P4 


2,  2,  2,  2 

Pi  -*2  ^3  ^4 


g 


exp 


i-»-  1 5/3  , -*■  1 5/3  i-»-  |5/3 

|P1-P2I  + I p3-p4 1 + lp’~  p- 


'1  4 1 


+ ip2-p3 


I 5/3  -*■  -*■  | 5 / 3 (-*■  -*■  | 5/3 

" I |J2_P4  " IP1-P31 


exP  ! Co  Ip-i~pJ)  + c c. ^ I P 2-P 4 1 ^ 


1 3 


(34) 


where  D is  the  — irradiance  diameter,  p"  is  the  i^  difference-vector  in 
g e 

the  transmitter  plane,  and  C is  the  log  amplitude  covariance  function. 

I 21 

It  may  be  shown  from  physical  reasoning,  or  from  tedious  but  straight- 
forward analytical  calculations,  that  the  covariance  functions  may  be 
dropped*  when  we  consider  the  limit  D^/p^^00. 


*It  may  also  be  shown  that  this  is  equivalent  to  dropping  the  amplitude 
term  in  the  perturbation  Green's  function,  as  done  in  Ref.  23.  In  the 
latter  effort,  the  phase  perturbation  was  assumed  to  be  gaussian,  and  the 
problem  under  discussion  here  was  reduced  to  a six-fold  integral  for 
numerical  evaluation.  Although  the  results  in  the  limit  of  large  Dg/P0 
agree  with  those  presented  here,  it  may  be  demonstrated  that  the  dropping 
of  the  amplitude  term  results  in  inaccuracies  at  smaller  values  of  Dg/pQ. 

D. 


We  let  x - p/D  , a = 
8 
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<*2> 


, and  write  (34)  as 


dx2dx^dx^  exp 
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23. 


V.  A.  Banakh,  et  al,  J.  Opt.  Soc.  Am.  64,  516-518,  April  1974. 
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We  now 

write 

lim  exp  | - (a)^3  f |x.-x  | 

a-x» 

1 

L 

+ l^3-x4|  5/3  + I x1-x4|  5/3  + |x2-x3!5/3 


5/3 


5/3 


y ^ I 5/3  1 

+ ixrx4'  J 


= 6/5ttF  (6/5)  (a)  2 exp  ^-(a)5/^3  |xj-x4 1 5^3  | 5 ( | x^-x^,  | ) 
+ 6/5tt T (6/5)  (a)  2expj  -(a)5^3  jx^-x^  <$  ( | x, -x,  | ) 


1 4 


i.e. 


(35) 


(36) 


6(|I  -i2l)-llm  5/6  ^rW5)  6,p  <c)5/3lVS2!5/3|  <37) 

a-Kn  ( ) 


where  we  note  that  this  satisfies  the  properties  of  a delta  function: 
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Using  (36)  with  (35),  we  have 
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Iterating  this  procedure  for  the  terra  in  each  of  the  integrands  of  the 
form  exp  j-(ct)5/3  l^l573!  =)6/5irr(6/5)a“2  6(1^-^!)  results  in 
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Similarly,  we  may  write 
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Using  the  delta  function,  this  becomes 

2 


lim 

0(-X» 


(0  0/3)  ^ j V"2 

■$)  (3/5) 

/kD\2 

1377  1 -if) 

Finally,  we  write  from  (AO)  and  (A2): 

2 


lim  a, 

a-*«° 


(*-)  ~ (3-)2  , 
<7«>2 


This  is  the  desired  result.* 


*In  order  to  make  the  development  of  the  present  section  agree  with 


Ref.  18,  it  is  necessary  to  multiply  I by  f*_ 
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that  of 


IV.  FUTURE  WORK 


In  a follow-on  effort,  we  will  obtain  more  complete  experimental 
results  on  short-term  statistical  effects,  using  the  analytical  understand- 
ings presented  in  Sec.  II,  and  we  will  carry  out  related  computer  simulation 

work. 

Comprehensive  experiments  on  finite-beam  target  illumination  effects 
will  be  completed,  with  varying  pathlengths  and  wavelengths.  In  addition, 
our  recent  analytical  advances  in  the  rigorous  determination  of  the  second 
moment  of  irradianc''  vs.  D/Pq  will  be  pursued,  with  supporting  numerical 

calculations . 
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